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X MM -Newton monitoring of X-ray variability in the quasar 

PKS 0558-504 * 
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Abstract. We present the temporal analysis of X-ray observations of the radio-loud Narrow-Line Seyfert 1 galaxy 
(NLS1) PKS 0558-504 obtained during the XMM-Newton Calibration and Performance Verification (Cal/PV) 
phase. The long term light curve is characterized by persistent variability with a clear tendency for the X-ray 
continuum to harden when the count rate increases. Another strong correlation on long time scales has been found 
between the variability in the hard band and the total flux. On shorter time scales the most relevant result is the 
presence of smooth modulations, with characteristic time of ~ 2 hours observed in each individual observation. 
The short term spectral variability turns out to be rather complex but can be described by a well defined pattern 
in the hardness ratio-count rate plane. 
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1. Introduction 

Active Galactic Nuclei (AGN) are variable in every ob- 
servable wave band. The X-ray flux exhibits variability on 
shorter time scales than any other energy band, indicating 
that the emission occurs in the innermost regions of the 
central engine. Therefore, a study of the X-ray variability 
is a powerful tool providing upper limits on the sizes of 
the emitting regions and the masses of central black holes 
and it allows one to probe the extreme physical processes 
operating in the inner parts of the accretion flow close 
to the event horizon. Although X-ray variability has been 
observed in AGN for more than two decades, its origin 
is still poorly understood. Narrow-Line Seyfert 1 galaxies 
(NLS1) often display rapid, large amplitude X-ray vari- 
ability as well as extreme long term changes (Forster & 
Halpern p96| , Boiler |ct al. 1997[ , Brandt |et~al. 1999] ), and 
therefore they represent the ideal class for an X-ray tem- 
poral analysis. PKS 0558-504 (z = 0.137, m B = 14.97) 
is one of the few radio-loud NLS1. In X-rays it is charac- 
terized by a steep spectrum [the photon index T ranges 
between 2.2 and 3.1, obtained from different instruments 
and energy bands (a detailed spectral analysis of XMM- 
Newton observations has been performed in a separate 
paper, O'Brien et al. 2001)], high luminosity [(2 — 5) x 



Table 1. EPIC PN individual observations. 



Observation date 


Duration 


Mode a 


Filter 6 


00/2/7 UT11:15-13:19 


7.2 ks 


FF 


M 


00/2/7 UT14:28-16:02 


5.4 ks 


LW 


M 


00/2/10-11 UT23:26-3:23 


14.0 ks 


FF 


T2 


00/3/2 UT18. 15-21:42 


12.4 ks 


FF 


M 


00/5/24 UT17:21-21:07 


13.6 ks 


SW 


M 



a Full Frame (FF), Large Window (LW), Small Window 
(SW). 

b Medium (M), Thin 2 (T2). 



10 45 erg s" 1 



with Hq 



- 1 Mpc" 



3o 



0.51 
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* Based on observations with XMM-Newton, an ESA Science 
Mission with instruments and contributions directly funded by 
ESA Member States and the USA (NASA). 



70 km s" 

and strong variability. A Ginga observation (Remillard 
[et al. 1991 ) showed an increase of the X-ray flux by 67% 
in 3 minutes, implying the presence of a relativistic flare. 
Recent ROS AT High Re solution Imager (HRI) observa- 
tions (Gliozzi |st al. 200C ) confirmed the presence of strong 
and persistent X-ray variability, which suggests the pres- 
ence of a rotating black hole, and ruled out external contri- 
butions to the high luminosity and variability from other 
nearby sources. 

In this paper we report the results of four XMM- 
Newton observations of PKS 0558-504 taken with the 
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European Photon Imaging Camera (EPIC) PN, which is 
the best suited instrument on board the ESA satellite for 
timing analysis purposes, due to its high time resolution. 
In Sect. 2 we describe the observations and data reduction. 
The long and short term X-ray variability are discussed in 
Sect. 3 and Sect. 4, respectively. Sect. 5 contains the main 
conclusions and a summary. 

2. Observations 

PKS 0558-504 was observed with the EPIC PN camera 
during orbits 30, 32, 42 and 84 (corresponding to February 
7 and 10, March 2 and May 24, respectively). The cam- 
era was operated in different observing modes and with 
various filters. In the Full Frame mode all the pixels of all 
CCDs are read out with time resolution of 70 ms and the 
full field of view (~ 30') is covered. In the Partial Window 
modes (Large Window and Small Window) a two dimen- 
sional imaging readout over part of the CCDs array is per- 
formed with increased time resolution (48 ms and 6 ms, re- 
spectively) . The primary reason for using filters is that the 
EPIC CCDs are not only sensitive to X-ray photons, but 
also to IR, visible and UV light. Therefore for objects with 
high optical to X-ray flux ratio, the X-ray signal might be 
contaminated by those photons. This is not the case of 
PKS 0558-504, for which different filters were used only 
for calibration purposes. Table 1 provides a summary of 
the individual observations. The data reduction and anal- 
ysis were performed using the last release of the XMM- 
Newton SAS (Science Analysis Software) package. From 
the available observations, the last reprocessed data sets 
without any indication of timing problems and only single 
events (i.e. events where the charge released by a photon 
is contained in a single CCD pixel) have been analyzed. 
The main reason for using only singles is related to the 
energy dependence in the relative fraction of "multiple" 
events. To quantify their contribution, a test with all valid 
events has been performed. The main effect is an increase 
of the count rate by ~ 10% and variations of the hardness 
ratio by ~ 20%, which would introduce additional uncer- 
tainties. Light curves were obtained by extracting pho- 
tons from circular region of radius of ~ 33" around the 
source center and subtracting the background taken from 
a source-free region (which amounted to ~ 1%). After con- 
firming the lack of significant variability on time scales less 
than minutes, the photons were binned into 200 s inter- 
vals providing a good signal to noise ratio. PKS 0558-504 
was chosen as a CAL/PV target with the purpose of in- 
vestigating the photon pile-up in the EPIC cameras. With 
this aim we compared the spectra from the central pixels 
with those from the outer parts of the PSF and found no 
significant indication for photon pile-up. 

3. Long term variability 

In Fig. we show the total light curve for PKS 0558- 
504 from February 7 to May 24 2000. Each data point is 
an average count rate over the individual observation in 
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Fig. 1. EPIC PN 0.2-10 keV light curve of PKS 0558-504 
for orbit 30, 32, 42 and 84. 
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Fig. 2. Hardness ratio (hard-soft)/(hard+soft) versus to- 
tal count rate. Every data point corresponds to 2000 s 
integration time. Soft band (0.2-1 keV), hard band (1-10 
keV). 



the range 0.2-10 keV. The indicated error of la, repre- 
sents the dispersion around the mean and gives an idea 
of the variability within each observation. The error bar 
is particularly large for the first data point because we 
merged two exposures (with very different mean count 
rates) taken during orbit 30 nearly one hour apart. The 
total light curve is characterized by a moderate variabil- 
ity, if compared with the ROSAT light curve where the 
count rate varied by a factor 5 in less than 3 days (Gliozzi 
|ct al. 200C| ). In the EPIC PN light curve the most promi- 
nent variation in the count rate is by a factor of 2 over a 
period of 20 days. 
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Fig. 3. Ratio of excess variances versus total mean count 
rate. Each filled circle is the average over the individual ob- 
servation. The triangles represent the average values over 
the two exposures during orbit 30. 



We sought evidence for spectral variability by comput- 
ing the hardness ratios as a function of the total 0.2-10 keV 
flux. As hardness ratio we define the difference between 
the count rates in the 1 keV < E < 10 keV band and 
those in the 0.2 keV < E < 1 keV band, divided by the 
0.2-10 keV count rate. These ranges where chosen to pro- 
vide a good signal to noise in both the soft and the hard 
light curves. The hardness ratios versus the total count 
rate with data from each observation binned in 2000 s in- 
tervals is shown in Fig. |[ The presence of a positive cor- 
relation was quantitatively tested by performing a linear 
least square fit, which confirmed the result at 5cr confi- 
dence level. On long time scales the spectral variability 
is correlated with flux variations such that the spectrum 
becomes harder when the count rate increases. 

Another interesting correlation was found between the 
variability in the hard band and the total count rate . The 
variability was studied using the excess variance (Nandra 
ct al. 1997), which is obtained by computing the variance 
of the overall light curve, subtracting the variance due 
to measurement error and dividing by the mean squared. 
Care must be taken in interpreting the excess variance of 
different observations as it is related to the length and the 



sampling of the time series (e.g. Leighly 1999a). However 
this drawback is circumvented in our analysis by consider- 
ing the ratio between the excess variance in the hard band 
and that in the soft range during each individual obser- 
vation. Fig. H shows a clear correlation between the hard 
to soft excess variance ratio and the mean count rate in 
the 0.2-10 band: the contribution of the hard band to the 
variability increases and becomes dominant as the mean 
count rate increases. From a direct inspection of the indi- 
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Fig. 4. Individual EPIC PN 0.2-10 keV light curves of 
PKS 0558-504 for orbit 30, 32, 42 and 84 with time bin- 
ning of 200 s. 

vidual values of the excess variance in the soft and hard 
band, this behaviour turns out to be due more to the ac- 
tual increase of the hard variability than to a depletion of 
the soft variability level. 



Based on ROSA T and ASCA data Fiore |et al. (1998) 
and Leighly ( 1999b| ) reported the presence of correlations 
between variability and X-ray luminosity and the steep- 
ness of the X-ray spectrum. A first analysis seems to indi- 
cate that the XMM data of PKS 0558-504 are consistent 
with the previous results, however a detailed study of such 
correlations is beyond the scope of this paper. A spectral 
study of the XMM observations of PKS 0558-504 can be 
found in O'Brien et al. (2001). 
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Table 2. Variability properties of the individual light 
curves. 



Orbit 


Count rate" 


b 

"rms 


soft 
v rms 


hard 
^rms 


30a 
30b 
30tot 


8.05 ±0.51 
12.65 ±0.78 
10.05 ± 2.39 


3.3 ± 1.2 
3.3 ±0.9 
55.0 ±3.0 


4.6 ± 1.8 

2.7 ± 1.0 
53.5 ±3.7 


0.4 ±2.0 
5.2 ±5.0 
61.0 ±14.6 


32 


10.99 ± 0.95 


7.0 ± 0.9 


7.9 ± 1.3 


7.1 ±4.0 


42 


16.30 ± 1.02 


3.6 ±0.5 


2.7 ±0.5 


9.7 ±3.0 


84 


10.54 ±0.51 


1.8 ±0.5 


2.0 ±0.7 


1.4 ±1.9 



a The errors on the count rate represent the dispersion around 
the mean. 



Excess variances are in 10 units. 



4. Short term variability 

The PKS 0558-504 EPIC PN light curves for each indi- 
vidual observations, binned into 200 s intervals, are shown 
in Fig. ||. Strong and persistent X-ray flux variability is 
observed during each observation. The most prominent 
increase in the count rate is observed during orbit 30 
(February 7), when the 0.2-10 keV flux increased by a 
factor 2 in less than 4 hours. On the other hand the most 
extreme temporal variation in the count rate occurred on 
February 10, when the count rate decayed by 18% in 1500 
seconds. This value is in agreement with that measured 
by the EPIC MOS during orbit 84 (O'Brien |st al. 2001] ), 
which led to the calculation of a lower bound for the ra- 
diative efficiency which exceeds the theoretical maximum 
for a non rotating black hole, under the assumption that 
the photon diffusion through a spherical mass of accreting 
matter is dominated by Thomson scattering. 

A visual examination of the individual light curves in- 
dicates that all of them show a similar long term vari- 
ability pattern with a common rise time of about 2 hours 
and, in some cases, small amplitude flares superimposed. 
The best method to quantify time variability without the 
problems encountered in the traditional Fourier analysis 
technique in case of unevenly sampled data is a struc- 
ture function analysis (e.g. Simonetti et al. 1985 ; Hughes 
|et al. 1992 ). The first-order structure function is the mean 
deviation for data points separated by a time lag r, 
SF(t) = ([F(t) - F(t + r)] 2 ). One of the most useful 
features of the structure function is its ability to discern 
the range of time scales that contribute to the variations 
in the data set: the characteristic time scales of the vari- 
ability are identified by the maxima and slope changes in 
the r — SF plane. For a stationary random process the 
structure function reaches a plateau state for lags longer 
than the longest correlation time scale. If a light curve con- 
tains cycles of period P, the SF will rise to a maximum 
at t = P/2 (Smith |et al. 1993| ). 
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Fig. 5. Structure functions of orbits 30, 32, 42 and 84. 
Time lags are in seconds. The dashed lines have been 
drawn to mark the interesting time lags interval. Most 
of the error bars are smaller than the symbol size. 

Fig. U shows that all the observations have either a rel- 
ative or absolute maximum around 2 hours. This means 
that a common typical time scale (which probably reflects 
the similar rise time in the light curves) characterizes all 
observations. In addition, during orbits 30 and 84 a fur- 
ther characteristic time scale around 3 hours seems to be 
present. A first confirmation for the presence of quasi- 
periodicity in the temporal behaviour of PKS 0558-504 
comes from a periodogram analysis, which yielded a strong 
signal at 2.3 hr. However, further longer observations are 
necessary for a firm confirmation. 

In order to seek evidence for spectral variability and 
for its origin on short time scales we have split each of the 
broad band light curves into a hard and a soft light curve 
and plotted the hardness ratios versus the time. While 
during orbits 32 and 84 (when the flux was at interme- 
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Fig. 6. Top panel: soft (0.2-1 keV) and hard (1-10 keV) 
light curves during orbit 30 ; data are binned in 200 s 
intervals. The soft light curve divided by the mean (empty 
squares) has been multiplied by 1.5, to avoid overlapping 
with the hard one (filled circles). Bottom: hardness ratio 
versus time with 200 s binning. 
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Fig. 7. Top panel: soft (0.2-1 keV) and hard (1-10 keV) 
light curves during orbit 42 ; data are binned in 200 s 
intervals. The soft light curve divided by the mean (empty 
squares) has been multiplied by 1.5, to avoid overlapping 
with the hard one (filled circles). Bottom: hardness ratio 
versus time with 200 s binning. 



diate values) no significant spectral variability was found, 
orbits 30 and 42 present an interesting and somewhat puz- 
zling spectral and temporal behaviour shown in Fig. ^ and 
Fig. [?|, respectively. The mean count rates and the excess 
variances in the soft, hard and broad bands are summa- 
rized in Table 2. 

During the first exposure of orbit 30 the variability, 
quantified by the excess variance ijj ms is strongly domi- 
nated by the soft photons (4.6 x 10~ 2 vs. 4.0 x 10~ 3 , for 
the soft and the hard band respectively), with the soft 
flux steadily increasing and the hard fluctuating around 
a constant value. As a consequence, while the total flux 
increases the hardness ratio decreases in time, with a 
slope of (—2.9 ± 0.7) x 10~ 2 /hr obtained from a linear 
least square fit. A totally different trend is observed dur- 
ing the second exposure, starting just one hour later. In 
this case the variability is dominated by the hard flux 
(er 2 ms : 5.2 x 10 -2 and 2.7 x 10~ 2 , for the hard and soft 
band, respectively), which increases faster than the soft 
emission and seems to peak before it. The hardness ra- 
tio is steadily increasing with the total flux and the time 
(slope: (1.2±0.4) x 10 _2 /hr). A similar spectral behaviour 
is observed during orbit 42: the variability is strongly dom- 
inated by the hard flux (<j 2 ms : 9.7 x 10~ 2 vs. 2.7 x 10~ 2 ) 
and the spectrum becomes harder when the total count 
rate increases (slope: (1.3 ±0.3) x 10 _2 /hr). 



On the basis of the five exposures available (two for 
orbit 30 and one each for orbits 32, 42 and 84), the spec- 
tral variability on short time scales of PKS 0558-504 can 
be summarized in the following way: when the average 
broad band count rate is at a low level (the mean value 
is around 8 cts s _1 ), the variability is dominated by the 
soft photons and the hardness ratio decreases as the to- 
tal flux raises. Conversely, when the average count rate is 
above a certain level (12.6 cts s _1 and 16.3 cts s _1 for the 
second exposure of orbit 30 and orbit 42, respectively), 
the variability is dominated by the hard band and the 
spectrum becomes harder when the count rate increases. 
At intermediate values of count rate (the mean is around 
10 cts s _1 both for orbit 32 and 84), the variability in the 
two bands is very similar and no significant spectral vari- 
ability is detected. 

5. Summary and conclusions 

We have presented the results of the EPIC PN observa- 
tions of the NLS1 PKS 0558-504 taken during the Cal/PV 
phase. We found: 

- During the XMM pointings PKS 0558-504 showed 
a moderate but persistent variability both on long 
(months) and short (minutes and hours) time scales. 
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— Strong correlations between hardness ratio and total 
flux have been found. When the mean flux is above 
a certain threshold, the spectrum becomes harder as 
the count rate increases. When the flux is below a cer- 
tain level, the spectral behaviour is in the opposite 
sense. At intermediate count rates no significant spec- 
tral variability is observed. 

— A further strong correlation between the broad band 
flux and the variability (expressed in terms of excess 
variance) in the hard energy band has been discovered: 
the contribution to the total variability of the hard 
band increases and becomes dominant as the mean 
count rate raises. 

— In all the observations the variability is characterized 
by smooth similar modulations, with typical time scale 
~ 2.2 hr, determined with the structure function anal- 
ysis. 

— The most extreme temporal variation of the count 
rate implies a radiative efficiency slightly larger 
than the theoretical maximum for accretion onto a 
Schwarzschild black hole. It is important to point out 
that this is a strictly lower limit, since during our anal- 
ysis only single events (which represent ~ 70% of the 
total X-ray photons collected) have been used. 

A detailed model of PKS 0558-504 which accounts for 
all timing and spectral features discovered is beyond the 
scope of this paper. Therefore, we limit ourselves to some 
general considerations. The most intriguing result is the 
tentative detection of a common characteristic time scale 
of ~ 2 hours. However a confirmation of this characteristic 
time scale has to await for longer observations, to sam- 
ple a higher number of cycles. A more detailed analysis 
to assess it s significance and robustness (e.g. as in Boiler 
et al. 2001 ) is deferred to a follow-up publication. This 
time scale might be associated to periodic phenomena oc- 
curring in the accretion disk, whose contribution seems 
to be dominant also on the basis of the spectral analy- 
sis of EPIC MOS data (O'Brien |et al. 200l[ ). Two possi- 
ble explanations which need to be quantitatively tested 
are: i) the presence of a hot spot orbiting the black hole 
and ii) the precession of the inner part of the accretion 
disk caused by misalignment of black hole and disk rota- 
tional axes. An alternative tempting hypothesis involves 
the presence of a jet: as PKS 0558-504 is a radio-loud ob- 
ject, beamed emission from a jet could partly contribute 
to the brightness and variability in X-rays. Radio-loud 
quasars are known to have flatter X-ray spectral indices 
than radio-quiet ones, and they are also known to be 



that also in one other of the few radio-loud NLS1, RX 
J0134. 2-4258, the spectral variability (characterized, as 
in our case, by a hardening of the spectrum during the 
flux increase) was tentatively interpreted as due to the jet 
emission (Grupe ;t al. 200"o| ). 



brighter X-ray sources (e.g. Yuan et al. 1998). A widely 
accepted explanation for this fact is that the X-ray emis- 
sion includes an additional harder component associated 
with the radio jet. In this scenario a naive interpretation of 
the temporal and the spectral variability associated with 
changes of the broad band flux, involves the presence of 
cither a precessing jet or moving knots with helical tra- 
jectories in a r elativ istic magnetized jet (Camenzind & 
Krockenberger 1992 ). At this point it is worth noticing 



The exciting results from the temporal analysis of a 
few short exposures demonstrate the extraordinary capa- 
bilities of XMM, in particular in light of the long obser- 
vations possible in virtue of the highly eccentric orbit. 
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